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OF  QUASI  STEADY  STATE  LOADINGS 
ON  THIN  ELASTIC  WINGS  IN  SUBSONIC  FLOW1 


by 

Sven  G.  Hedman 


SUMMARY 

Wing  systems  that  contain  several  wing  surfaces  are  considered.  The  surfaces  arc  arbitrarily  divided 
in  the  chordwise  and  spanwisc  directions  into  panels  (“boxes”).  The  panel  load  is  simulated  by  a 
horseshoe  vortex,  and  the  boundary  condition  is  fulfilled  on  every  panel  at  one  point.  For  each  surface 
the  vortices  and  the  control  points  arc  positioned  on  a  wing  chord  plane.  These  planes  are  perpendicular 
to  the  ij-  or  r-axis.  The  incidence  distribution  or  the  load  distribution  may  be  found  from  a  system 
of  linear  equations,  when  one  of  these  distributions  is  given.  The  calculations  can  be  performed  for 
clastic  wings,  if  the  clastic  properties  are  known  in  the  form  of  a  deformation  matrix. 

The  calculation  procedure  has  been  programmed  for  an  electronic  computor.  Good  correlation  has 
been  obtained  between  the  results  from  this  method  and  those  from  other  lifting  surface  methods  or 
experiments. 


Stockholm,  October  1965 


1  The  work  was  initiated  by  the  author  in  1900,  while  he  was  employed  by  the  Boeing  Co.  It  lias  been  ex¬ 
tended  at  the  FFA  and  partly  sponsored  by  the  SAAB  Co.  A  similar  method  has  been  developed  independently  by 
Rubbcrt,  Ref.  1. 
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Z  U  SAM  M  E  NFAS  SUNG 

Es  wcrdcn  Tragflugclsystcmc  bctrachtet,  die  aus  mchreren  Flugclflachcn  bcstehen.  Die  Flaclicn 
wcrden  in  Profilsclincnrichtung  und  in  Richtung  der  Spannwcite  bclicbig  in  Felder  (boxes)  unter'.ilt. 
Die  Fcldbclastung  wird  durcli  eincn  Ilufciscnwirbel  nachgebildct  und  die  Randbedingungcn  wcrden 
in  jedcm  Feld  an  cincm  Punkt  crfulit.  An  jcder  Flache  sind  die  Wirbcl  und  die  Aufpunkten  auf  cincr 
Flugelschncnflaclic  angcordnet.  Dicse  Flaclien  stchcn  senkrcclit  zur  odcr  z-Aclisc.  Die  Anstcll- 
winkclvcrteilung  oder  die  Lastverteilung  kann  aus  eincm  System  linearer  Glcichungcn  gcfundcn  wcr¬ 
den,  wenn  cine  dicscr  Vcrtcilungcn  gegeben  ist.  Die  Bercchnungen  konnen  fur  clastisclic  Flugcl  durch- 
gefuhrt  wcrden,  sofern  die  elastischcn  Eigcnscliaftcn  in  Form  eincr  Formandcrungsmatrizc  bekannt 
sind. 

Der  Reclienvorgang  wurdc  fur  cinen  Rechenautomat  programmiert.  Die  Ergcbnissc  dicscr  Mc- 
thode  stimmen  gut  mit  jenen  andercr  Tragflachenbcrechnungen  odcr  Vcrsuchen  ubcrcin. 


MfiTHODE  A  LATTIS  DE  TOURBILLON  POUR  LE  CALCUL  DE  CHARCES 
QUASI  STABLES  SUR  AILES  DEFORMABLES  ET  MINCES 
EIN  REGIME  SUBSONIQUE 


par 

Sven  G.  Hedman 


RESUME 

Lcs  systdmes  de  voilure  comportant  plusieurs  surfaces  d’ailcs  sont  pris  ici  cn  considdration.  Les 
surfaces  sont  arbitrairement  divisdes  longitudinaleincnt  ct  transvcrsalemcnt  en  panneaux.  La  charge  du 
panneau  cst  simuldc  par  un  tourbillon  lie  ct  deux  tourbillons  libres  ct  la  condition  liinitc  cst  appliqude 
en  un  point  de  cliaque  panneau.  Pour  cliaque  surface,  les  tourbillons  ct  les  points  de  contrdle  sont 
positionnds  scion  un  plan  rdfdrence  de  1’aile.  Ccs  plans  sont  perpendiculaircs  aux  axes  y  on  z.  La  di¬ 
stribution  d’incidence  ou  la  distribution  dc  charge  peuvent  Gtre  trouvdes  grdee  A  un  systGme  d’dquations 
lindaires  fi  condition  quo  l’une  des  deux  distributions  soit  donnde.  Les  calc'.ls  peuvent  etre  cffcctues 
pour  des  ailcs  ddformables  si  les  propridtds  d’dlasticitd  sont  connucs  sous  forme  d’une  matricc  dc 
ddforination. 

Les  proeddds  de  calcul  ont  dtd  programmes  pour  unc  calculatricc  dlectroniquc.  Unc  bonne  corrdlation 
existe  entre  lcs  rdsultats  obtenus  par  ccttc  indthode  ct  ceux  ootenus  par  d’autres  mdthodcs  ou  des 
essais. 
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1.  INTRODUCTION 


In  the  literature  a  great  number  of  meth¬ 
ods  arc  available  for  the  approximate 


in  subsonic  flow.  Most  of  these  methods 
were  constructed  before  the  arrival  of  the 


automatic  coinnutor,  and  thus  the  number 
of  unknowns  had  to  he  kept  low.  This  was 
accomplished  by  assumptions  on  the  chord- 
wise  and  spanwisc  load  distributions.  These 
assumptions  were  well  founded  for  a 
majority  of  wing  designs.  One  of  these 
methods  was  the  vortex  lattice  method  of 


Faulkner.  However,  with  the  electronic 
computor  available,  the  number  of  un¬ 
knowns  in  a  system  of  equations  can  be 
permitted  to  rise  considerably,  and  therefore 
the  assumptions  on  load  variations  arc 
unnecessary,  and,  as  they  ore  not  always 
suitable,  it  is  belter  to  avoid  them. 

This  report  describes  an  elementary 
method  for  load  calculations.  Let  the 


continuous  vortex  distributions,  which  re¬ 
present  the  continuous  loading  on  the  lifting 
surfaces,  be  replaced  by  a  system  or  discrete 
vortices.  The  strengths  of  the  vortices  arc 
determined  bv  the  requirement  of  tangen¬ 
tial  (low  in  as  many  points,  control  points, 
as  there  arc  vortices.  The  system  can  be 
thought  of  as  a  collection  of  horseshoe 
vortices.  Each  horseshoe  vortex  induces  the 


same  flow  field  as  an  elemental  lifting  area 
of  the  wing,  a  panel,  does.  For  each  wing 
surface  the  vortices  and  the  control  points 
arc  positioned  on  a  wing  chord  plane.' 

Because  the  panels  can  be  laid  out  in 
rather  arbitrary  patterns,  the  boundary  con¬ 
dition  can  be  fulfilled  on. any  planfornl. 


However,  the  proposition  that  the  vortices 
arc  arranged  in  planes,  limits  the  method 
to  problems  with  such  incidences  that  no 
vortex  roll  up  exists  over  the  leading  edges, 
and  that  the  trailing  vortex  sheets  arc  close 
to  their  original  wing  chord  planes. 


2.  NOTATIONS 

b  {/-coordinate  of  the  leading  edge  of  whig 
part,  26max  equals  span 
CD  drag  coefficient 
CL  lift  coefficient 
C,  rolling  moment  coefficient 
Cm  pitching  moment  coefficient 
Cp  local  lift  coefficient 
C„  yawing  moment  coefficient 
Cy  side  force  coefficient 
c  wing  chord 

h  half  the  width  of  horseshoe  vortex 
spanwisc  loading  coefficient 
P,  Q,  R  aerodynamic  influence  coefficients 
p  rate  of  roll 

q  dynamic  pressure  or  rate  of  pilch 
s  area  of  wing  panel 

S,  matrix  of  aerodynamic  influence  coefficients 
u,  v,  w  components  of  disturbance  velocity 
along  x,p,z-ax\s,  respectively 
V  free-stream  velocity 
x,  ij,  z  Cartesian  coordinates 
a  wing  angle  of  attack 


z.w 


c 
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P  l  l  -M-,  where  M  is  Mach  number 

r  vortex  intensity 

<5  wing  camber 

0  2q-(xD-xe)/c„t 

7.  angle  of  sweep  of  bound  vortex 

o  density  of  air 

^  PUl^CJLX 

a)  a  +6  +0  panel  angle  of  attack 

SUBSCRIPTS 

D  control  point 
V  vortex 

3.  FLOW  MODEL 

The  lifting  system  may  consist  of  several 
wings,  or  it  can  be  the  horizontal  projection 
of  a  complete  airplane. 

For  the  derivation  of  the  formulas  relating 
incidence  and  load,  the  lifting  surfaces  will 
he  replaced  by  a  number  of  plane  quadri¬ 
laterals  or  panels.  Thus  the  planform  of 
the  original  system  will  he  approximated 
hv  polygons,  and  the  angle  of  attack  distri¬ 
butions  of  the  surfaces  will  be  expressed 
through  the  incidences  of  the  plane  panels. 

The  load  carried  by  one  of  these  panels 
induces  a  flow  field  that  can  be  calculated 
with  the  aid  of  a  horseshoe  vortex.  The 
bound  vortex  will  he  positioned  on  the 
panel’s  quarter  chord  line  for  reasons  out¬ 
lined  in  Appendix  1.  The  control  point, 
where  the  condition  of  tangential  flow  is 
satisfied,  will  he  positioned  on  the  three 
quarter  chord  line  of  the  panel  and  halfway 
between  that  line’s  inboard  end  and  out¬ 
board  end.  For  each  wing  the  vortices  and 
control  points  will  he  in  a  representative 
wing  chord  plane,  perpendicular  to  the 
y-axis  for  a  vertical  tail  and  to  the  :-axis 
for  a  wing.  The  free  vortices  are  taken 
parallel  to  the  x-axis. 

At  every  control  point,  the  normal  velocity 
induced  by  all  the  vortices  must  cancel  the 
normal  component  of  the  free  stream  flow. 
The  induced  velocities  are  calculated  ac¬ 
cording  to  Biot-Savart’s  law’.  As  the  free 
vortices  are  assumed  to  he  parallel  with 
the  x-axis,  the  equations  for  the  induced 
velocities  will  he  linear  in  the  vortex  inten¬ 
sities. 


4.  GEOMETRY  OF  THE  WING 
PANELS 

The  method  is  not  very  sensitive  to  the 
pattern  of  panels  chosen  to  represent  the 
wing,  except  that  panels  behind  each  other 
should  he  in  streamwisc  columns.  Otherwise 
a  control  point  of  one  panel  may  lie  very 
close  to  the  trailing  vortex  of  another  panel, 
where  the  induced  velocity  is  high  and  not 
representative  of  the  average  induced  ve¬ 
locities  in  the  range  between  the  trailing 
vortices. 

The  horseshoe  vortex  is  geometrically 
defined  by  the  coordinates  x>,  ijv,  zv,  h,  7, 
and  a  control  point  by  xD,  yp,  zD,  as  shown 
in  Fig.  1. 


Fig.  1.  Wing  divided  into  panels  and  panel  with  vortex. 


Rather  than  determining  all  these  coor¬ 
dinates  manually,  it  is  advisable  to  define 
the  input  for  the  calculations  through  the 
major  dimensions  of  the  wing  and  the 
number  of  columns  and  rows  of  panels. 
The  scheme  in  Appendix  2  may  then  be 
follow'ed  for  the  determination  of  the 
coordinates  for  vortices  and  control  points. 

The  slope  of  a  panel,  co,  is  composed  of 
several  items.  The  stationary  contributions 
arc  the  reference  angle  of  attack,  a,  which 
is  assumed  to  have  one  value  for  the  whole 
lifting  system,  and  the  camber  distribution 
or  aileron  deflection  6. 

Quasi  steady  contributions  appear  for  a 
rotating  airplane.  Let  0  be  the  angular 
velocity  of  rotation  around  the  axis  x  =  x0, 
z  =  0.  At  a  point  xc,  ijD,  zD  the  normal  ve¬ 
locity  due  to  rotation  is  ~0-(xD  ~xfl).  The 
corresponding  load  is  the  same  as  for  a 
stationary  surface  so  cambered  that  the 
normal  velocity  is  -0-(xc-xo).  Let  0  be 
the  panel  incidence  that  corresponds  to  a 
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rotation  around  the  p:.tch  axis  (in  x  =  xfl), 
and  <f>  the  incidence  corresponding  to  a 
steady  roll. 


through  summation  over  all  the  vortices. 
Let  there  be  n  of  them  on  planes  parallel 
to  the  x,  y-planc. 


0-(xD-xe)  xD-xBn 
V  q  c„f  “ 


0) 


where  6  and  ij,  are  positive  for  a  sinking 
trailing  edge  and  a  sinking  starboard  tip 
respectively.  <7  and  p  are  the  dimensionless 
rates  of  pitch  and  roll  resp.,  and  fr„.T  is 
the  wing  half  span. 

f=fiW(210 

P  =  4-bvzxx!V 


The  total  panel  incidence  to  is  obtained 
through  summation 


a)  —  a-rS~0  —  <}> 


Angles  are  positive,  when  dzjdx  <  0. 


5.  AERODYNAMICS  OF  A  SET 
OF  PANELS 

The  expression  for  the  vortex  induced 
flow  field  is  derived  with  the  aid  of  Biot- 
Savart’s  law  which  implies  direct  applica¬ 
tion  only  to  incompressible  flow.  However, 
if  the  Prandtl-Glauert  rule  is  accepted,  the 
validity  of  the  derivation  can  be  extended 
to  all  subcritica!  Mach  numbers. 

Appendix  3  is  a  derivation  of  the  for¬ 
mulas  which  give  the  induced  velocity 
components  u,  v,  w  in  incompressible  flow 
at  a  point  (xD,  f}tjD,  fizD)  due  to  a  horseshoe 
vortex  strength  F;  at  (xv,  fiijv,  ftzv). 

Vu-Qv  r,  (2) 

wu~itvr, 

The  influence  coefficients  P,  Q,  R  depend  on 
the  geometry  of  the  horseshoe  vortex  and 
on  the  position  of  the  downwash  point 
relative  to  the  vortex.  The  coefficients  are 
calculated  with  Eq.  (A  3.1).  The  total  in¬ 
duced  velocity  components  are  obtained 


«i=  2  uu=  2  Pyr> 

1-1  j-i 

•  «.=  2  »„=  2  Qt,r,  (3) 

i-i  1-1 

w,=  2  2W 

1- 1 

A  cruciform  configuration  may  carry  load 
on  both  horizontal  arid  vertical  surfaces. 
Let  m  —  n  be  the  number  of  vertical  surface 
vortices.  To  avoid  introducing  a  further 
set  of  equations  to  calculate  the  influence 
coefficients  for  these  vortices,  an  x'.y'.z'- 
coordinale  system  will  be  inserted  accord¬ 
ing  to  the  sketch  below. 


Uy=uw  =  Pi  1  r  j 

v.f  =  —  wu  =  —  R,i  r t 
;o,j  =  vtj  —  QnVj 


11  +  1  <  /  <  m 


The  primed  influence  coefficients  are  ob¬ 
tained  from  Eq.  (A3.1),  if  y  and  z  from 
the  configuration  are  inserted  into  the 
equations  in  the  places  of  —z  and  y  respec¬ 
tively.  The  total  induced  velocity  for  a 
cruciform  configuration  is  obtained  through 
summation  over  all  the  m  vortices. 


«.=  I  F„r,+  2 

I- 1  f-n+l 

v,=  2  Q„T,-  2  K,T, 

I- 1  >-n+l 

W.-ZW+  I  Ql,r, 

/-I  I~n\ 1 


s 
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The  condition  of  tangential  flow  at  the 
control  points  is  fornuilalcd  under  the 
assumption  that  the  surface  slope  fho,  is 
small. 

—fat  =■%'  -  y  i<i/j  +  i  «'<« 

I  f  J-l  1  I  /- n*l  *  * 

II  B  r  n  ri 

-M=Vf=20(/y-  2  ^r/, 

*l  j-i  *:  j-n»i  >i 

n+l=Si</n 

T  will  he  replaced  by  CB,  where  Cp  is 
the  coefficient  for  the  average  pressure 
difference  between  the  panel’s  lower  and 
upper  surfaces.  The  arca.ofJhe  transformed 
-panel  is  '2 fill  -2(.td  -.rr). 

C  fincomnrl  ~ _ gl  Ft'  -fas 

Cp^incompr)  - 

=E / _ I _ Yi 

V  XD,-Xy,  V 

The  compressible  flow  pressure  coefficient 
CP  is  obtained  from  Prandtl-Glaucrt’s  rule. 

r  -1  rJ  1  V> 

The  tangential  flow  equations  arc  then 
written  in  terms  of  pressure  coefficients 

V  f  n  V 

-  co  i*=Pyt  |^2  l*i ,(xD,  -  Xri)  Cn 

+  2  Qll(Xffj—Xy/’)  yr  Cp/j  ,  l</l 
i-Hl  *J  J 

V  [  "  V 

— co  i =p  Yt  ^2  0  h(.xD) — Xv,)  Yt  Cpi 

-  2  K  (xdi  -  xvi)  X  c J ,  n  + 1  <  I  <  /n 

J-n+1  J 

In  a  parallel  flow  field  of  constant  velocity, 
V,  =  Vj  =  V.  For  this  condition  these  equa¬ 
tions  can  be  reduced  to  a  matrix  equation. 


{a)}-— ^[SJ  .T0~av  {Cp};  [SJ- 


The  upper  parts  of  the  coiumns  of  co  and 
Cp  contain  values  for  the  horizontal  surfaces 


and  the  lower  parts  values  for  the  vertical 
ones.  With  the  above  equation,  the  camber 
to  can  be  calculated  for  a  wing  that  carries 
the  loading  Cp  in  a  compressible  subsonic 
flow. 

However,  for  an  clastic  wing,  the  camber 
will  change  when  the  wing  is  unloaded, 
'flic  clTccls  of  elasticity  are  treated  in 
Appendix  4. 

The  load  distribution  that  corresponds  to 
a  given  incidence  distribution  can  be  cal¬ 
culated  after  inversion  of  Iiq.  (4). 


{«}  (5) 


After  the  pressure  distribution  is  known, 
force  and  moment  coefficients  are  easily 
summed. 

n 

2  s:  C' pi 

(«) 


2  S1  ’  CXrtf  x\’i)  Cp/ 

,)-l _ 

n 

Crct  2  SI 

l-l 


2  s)Cp/ 


/•  1+1 
W - iT 


2  S1 '  C^rcf  ~  *V,)  Cp/ 

P  l-n  +  1  _  _ 

C-'n  n 

2  Sj 

l-l 

n  m 

2  sjiJ Vj  Cp/  +2  s; ~Vj  Cp/ 

C,=>sl - tS*l (8) 

2tm«  2  Sj 

J-l 

In  these  equations,  s  is  the  panel  area 

s,=41il(xDl-Xyl),  (9) 

xrc(  is  the  reference  point  for  pitching  and 
yawing  moments,  cI0l  is  the  reference  length 
for  the  pitching  moment  coefficient,  and 
2 hm„  (“total  span)  is  the  reference  length 
for  the  i  filing  moment  coefficient. 
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For  many  problems,  the  spamvisc  load¬ 
ing  coefficient  for  unit  lift  coefficient  is  of 
great  interest.  This  parameter  is  commonly 
written  (<Vc)/(CL-ciy),  where  c,,  is  the 
average  chord.  It  can  also  he  defined  as  the 
relative  change  in  lift  due  to  a  change  in 
the  dimensionless  spamvisc  coordinate. 

.,  A(Llg) 

'  L'.q  'Ay 

t 

bmxx  2  (.XD/  —  *Ij)  Cpj 
l-l  — ~ - 

2  ■  (• xVf)  Cpj 

1-1 

where  s  <  /  is  the  interval  in  j  corre¬ 
sponding  to  the  panels  in  column  /. 

From  the  spamvisc  loading,  the  drag 
coefficient  is  readily  obtainable.  Ref.  2 
contains  a  suitable  expression  for  symmet¬ 
rically  loaded  wings. 


well  vortex  as  control  point  are  situated  on 
the  starboard  side.  The  induced  velocities 
arc  written  analogously  to  Fq.  (2) 

!.r=irir 

<»=/$•»  rj» 

The  contribution  of  the  port  side  vortex  to 
the  velocities  at  the  starboard  control  point 
are  the  same  as  the  starboard  side  vortex 
velocity  contributions  in  a  control  point 
placed  symmetrically  on  the  port  wing 
(except  for  the  sign  of  the  y-componcnl). 
The  latter  approach  is  chosen  because  it  is 
simpler  to  move  the  control  point  than  the 
vortex  from  starboard  to  port  wing  half. 

"ir’ = W  r?> = p\fn  rj“ 

= q«p>  =  -  Q<f‘>  IT 

«$» = ff{;»  rjw = lift1'  I?’ 


it  2  si 

CvlCi = •  [Lf + U?  v  L? + 


-AJ- (0.0561  LJ-r  0.7387  L0 
-  Li  •  (0.7352  L;  +  0.8445  L,)] 


L[,  /„«,  LJ,  and  LJ  arc  the  values  of  L'  at  the 
spamvisc  stations  y,/hmtx  =  0.9239,  y2/hmiX  = 
0.7071,  y3/hm„  =  0.3827,  and  y4/hmxx  =  0,  re¬ 
spectively. 

So  far  no  restrictions  or  symmetry  have 
been  made.  If  the  problem  is  symmetric 
or  antisymmetric,  great  savings  in  the 
amount  of  computation  for  a  given  accuracy 
can  be  made.  For  planrorms  that  are 
symmetric  with  respect  to  the  .c,  z-planc, 
only  one  half,  say  the  starboard  half,  of 
the  configuration  needs  to  be  defined  geo¬ 
metrically. 

When  the  camber  distribution  is  sym¬ 
metric,  the  vortex  intensity  rjs)  of  the  y-lh 
panel  of  the  starboard  side  is  the  same  as 
the  vortex  intensity  r,(p>  of  an  equally  posi¬ 
tioned  pane!  of  the  port  side.  The  induced 
velocity  due  to  such  a  panel  pair  is  the 
sum  of  each  panel’s  contribution. 

Let  the  superscript  (ss)  indicate  that  as 


The  influence  coefficients  with  superscript 
(ss)  arc  obtained  from  Eq.  (A3.1),  and  the 
ones  with  superscript  (/»)  can  also  be 
obtained  from  that  equation,  if  ijD  in  the 
equation  is  replaced  by  (— yD). 

The  total  velocity  components  arc  arrived 
at  by  summation. 

2  (JT + Pip)  r, 

i-i 

"i-icor-ODr, 

l-l 

wi-iw+TC 

l-l 


n  is  the  number  of  panels  on  one  wing 
half.  The  Eqs.  (4-7)  can  be  used  to  solve 
symmetric  problems,  if  the  following  state¬ 
ment  is  met: 

p^pir+pir 

Q^Qir-oir 

p^pir+nir 

Let  finally  the  camber  distribution  be 
antisymmetric,  and  apply  the  same  reas¬ 
oning  as  before.  It  should  be  recognized 

that  Fjp)  =-  r;j). 


10 


FFA  REPORT  105 


U-l^-W-rj 

I  " 

*»,=  5  (or-oro-r, 

i-i 

«««=  i  US1*- *«■*)- rf 

l  M 

Antisymmetric  cases  can  then  be  solved 
by  Eqs.(4,  5)  and  (S).’ 

r  Pu-Pff>-Fg» 

]  Qv-QF+Qg* 

\  nls=it?>-Jsr> 


6.  EXAMPLES  OF  CALCULA¬ 
TIONS 

The  presented  method  has  been  program¬ 
med  for  SAAB’s  computor  D21.  A  large 
number  of  calculations  has  been  performed. 
Some  of  them  which  arc  of  general  interest 
will  be  discussed. 

Fig.  2  contains  a  comparison  of  the 
spamvisc  distributions  of  load  and  center 
of  pressure  calculated  by  the  lifting  line 
method,  by  several  lifting  surface  methods, 
and  measured  by  Malavard  in  an  electric 
potential  tank.  The  wing  is  a  45°  delta  wing 
of  aspect  ratio  3.  The  SAAB  calculations 
were  performed  on  a  model  with  four  rows 
of  panels  distributed  over  the  chord  and 
with  10  columns  on  the  half  span.  Material 
for  the  comparison  is  taken  from  Refs.  3 
and  4.  The  table  in  the  figure  shows  the 
integrated  values,  lift  coefficient  and  total 
chordwise  center  of  pressure.  The  agree¬ 
ment  is  good. 

To  investigate  how  much  the  panel  pat¬ 
tern  influences  the  calculated  values  of  lift 
coefficient  and  center  of  pressure  several 
patterns  have  been  chosen  as  input  for  this 
particular  delta  wing.  In  all  eases  the  chords 
were  divided  into  equal  rows  of  panels, 
and  in  the  five  first  eases  the  half  span  was 
divided  into  equally  wide  columns  of  pan- 


fc - 


Lifting  line  3.6S  0.532 

Lifting  surface  Vcissinger  —  -  —  3.01  0.521 

Lifting  surface  Molthopp  -  3.06  0.535 

Lifting  surface  Falfcner  3.19  0532 

Potential  tank  Malavard  3.10  0535 

Vortex  lattice,  present  method  O  3.14  0531 

Fig.  2.  Loading  for  a  delta  wing.  Comparison  between 
methods. 

els.  In  the  sixth  ease  the  spanwisc  division 
was  made  in  intervals  of  the  angle  ar/20, 
ij  >=  sin  (nn/20),  0  <  n  <  10. 


Case 

Number 
ot  rows 

Number 
of  columns 

Lift  curve- 
slope 

Center 
ol  pressure 

x-coord 

«/-coord 

1 

1 

5 

3.20 

0531 

0.437 

O 

2 

5 

3.20 

0530 

0.430 

3 

4 

5 

3.20 

0538 

0.43G 

4 

S 

5 

3.20 

0539 

0.435 

5 

4 

10 

3.14 

0534 

0.429 

G 

4 

10 

3.17 

0534 

0.428 

All  values  arc  within  2  %  of  one  another. 


Measured  and  calculated  loadings  arc 
•mpared  in  Fig.  3  for  a  delta  wing  con¬ 
figuration  in  a  M  =  0.7  air  stream.  In  the 
calculations  the  body  was  not  considered. 
Instead  the  exterior  wing  was  extended  to 


DOWNWASH  AIIW.E  C1 
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Fig.  3.  Measured  and  calculated  loading  on  a  delta  wing  at  M  -  0.7. 


Fig.  1.  Measured  ami  calculated  conditions  away  from 
the  surface  of  a  swept  wing  at  8.2°  angle  of  attack. 


the  plane  of  symmetry.  The  loading  al 
i]  -  0.17  is  overestimated  by  the  theory. 
For  the  rest  of  the  wing  the  agreement 
between  measurements  and  calculations  is 
belter.  However,  this  agreement  deteriorates 
with  increasing  spanwisc  coordinate  and 
incidence.  This  discrepancy  is  due  to  the 
increasing  strength  of  the  free  vortex  over 
the  leading  edge. 

The  flow  conditions  outside  the  wing 
chord  planes  can  also  be  determined.  The 
calculated  data  are  shown  together  with 
experimental  ones  in  Fig.  4.  From  the 
experiments  reported  by  Alford  in  Ref.  5, 
the  effect  of  S.2°  angle  of  attack  was  ob¬ 
tained  by  subtracting  values  for  zero  angle 
of  attack  from  values  for  S.2°.  The  com¬ 
parison  is  limited  to  one  line  ij  =  3b/4; 
z=-0.15c,  where  b  is  the  wing  semispan, 
and  c  is  the  local  chord.  In  points  along 
this  line  the  disturbance  velocities  u,  v  and 
w  were  calculated.  They  are  shown  in  the 
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sr 

»«  3JCD*  t«E,  BX.Vy.5V, 


PI  A—* 

Fiji  5.  Measured  and  calculated  lift  coefhricJiis'fw  a 
wing  tu  a  propeller  slipstnrun  at  15*  angle  of  attack. 

graphs  as  a  local  dynamic  pressure  qjqa. 
sidewash  angle  a,  and  dowmvash  angle  e. 
All  the  flow  characlcrislics  arc  well  pre¬ 
dicted  by  the  calculations  except,  curiously , 
the  upwash  angle  under  the  wing. 

Finally.  »nc  example  with  a  wing  in  the 
slip  stream  of  a  propeller  will  lie  shown. 
The  experimental  data  arc  taken  from  Ref. 
G.  The  approximate  layout  of  the  wing  and 
propeller  appears  in  Fig.  5.  The  panel 
pattern  is  chosen  such  that  some  panels 
arc  completely  contained  inside  the  slip¬ 
stream  at  a  thrust  coefficient  C'T  =  6.14  and 
the  other  pancis  arc  exterior  to  the  slip¬ 
stream.  The  wing  and  propeller  axis  form 
an  angle  a  to  the  free  stream  direction,  a 
flap  covering  all  the  trailing  edge  is  de¬ 
flected  an  angle  6.  The  boundary  condition 
is  fulfilled  in  a  plane  parallel  to  the  free 
stream  velocity  V.  VP  is  the  propeller  slip¬ 
stream  velocity.  For  a  panel  situated  in  the 
slip  stream  the  propeller’s  dowmvash  is 
considered,  and  only  the  vector  w  is  con¬ 
tributed  by  the  vortex  system.  For  small 
angles  iv  becomes 

iv  —  (cc-r  £)(  V  4-  Fp)  —  aVp 

Fig.  5  also  contains  a  chart  that  shows 
Cl(6,  Ct)  for  a  =  15°,  C'T  is  the  propeller 


thrust  coefficient  based  on  free  stream  ve¬ 
locity  and  wing  area.  It  appears  that  a  good 
estimate*  sFUie  lift  eoc-fficient'can-bc  made, 
when  no  flow  separations  arc  present.  Flap 
separations  should  increase  with-  flap-  de¬ 
flection  and  decrease  with  thrust  coefficient. 
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APPENDIX  1 

LOCATION  OF  VORTICES  AND 
CONTROL  POINTS 

In  the  two-dimensional  ease  the  vortices 
and  control  points  can  be  so  located  on  a 
straight  airfoil  chord  that  the  vortex  lattice 
gives  the  same  lift  and  pitching  moment  as 
thin  airfoil  theory  does. 

Let  the  chord  be  divided  into  r  equal 
parts,  and  let  the  vortex  be  placed  a  distance 
£  times  the  part  chord  from  its  leading  edge 
and  the  control  point  a  distance  q  times 
the  part  chord  behind  the  vortex.  The 
following  figure  applies  to  the  i-th  part. 
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The  x-coordinate  for  Ihc  i-lh  control  point 
is  xft. 


After  identification  of  the  two  equations 
for  L  and  M  respectively,  y  and  -  can  he 
solved. 


The  inversions  have  only  been  performed 
for  r  =  2,3. 


The  j'-th  vortex  is  positioned  in  xVf. 


i)+f] 


The  velocity  induced  at  xD.duclo  the  vortex 
at  xVt  is  urtJ. 

r ,-r 

W,)  2.T c(«-i+,) 

The  distribution  of  induced  velocities  is 
related  to  the  vortex  distribution  by  a  set 
of  linear  equations. 

w-srlRTdm 

Insert  this  equation  into  the  equation  for  the 
total  lift  of  the  airfoii. 


r=2;  L-1-  .WJ'-r-ri,*1)] 

L«-y-rj?j  ;  Uo?-  *)• 1  -in 
=  3;  I  "1 

1 )  y(rf-  2) 

4  *»  2  ’ 

4 

w_n. 


r,(y-i)-(y-2)  qOr-QOj--) 
4  ’  2  ’ 

Of-I)(if~0 

4 


L=„v-iij{r}=5rH25.ii][r-~]"l{u;} 

Likewise  an  expression  for  the  moment  is 
obtained 

JZ-gK-  lav]{n=olr~-iy-  14-fj 


From  thin  airfoil  theory  the  lift  and 
moment  arc  known. 


L  =  In  V’2  •  c  •  2rr  •  wj  F 
.1/  =  L’c/4 


The  solutions  for  and  £  arc 

y  =  ■?  =  i 

if  r  =  1,  2,  3.  These  y-  and  f- values  numer¬ 
ically  satisfy  Fq.  (A  1.1)  for  useful  values 
of  r. 

The  same  vortex  and  control  point  posi¬ 
tions,  i.c.  |  =  },  y  =  1,  have  been  chosen 
also  for  the  three-dimensional  ease.  Calcu¬ 
lations  with  different  panel  patterns  for  one 
wing  planform  have  given  very  consistent 
values  of  lift  coefficient  and  center  of  pres¬ 
sure,  which  is  shown  in  the  table  on  P.  10 
and  in  Fig.  2.  The  suggested  locations  arc 
therefore  thought  to  be  very  useful  approxi¬ 
mations. 
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APPENDIX  2 

DIVISION  OF  THE  LIFTING 
SURFACE  INTO  PANELS 

Let  the  real  planform  be  approximated 
with  a  number  of  trapezoids.  In  Fig.  A2.1 
three  trapezoidal  wing  parts  represent  star¬ 
board  half  of  the  wing.  The  parts  arc 
defined  by  the  coordinates  a,  b,  c.  Each 
wing  pari  is  divided  in  the  spanwisc  and 
chordwise  direction  to  form  q  columns  and 
r  rows  of  panels. 

Fig.  A2.1  shows  that  the  interval  bk  <  y  < 
ft**,  is  divided  at  stations  y  =  ckt,  and  that 
the  chordwise  division  is  made  by  straight 
lines  between  the  points  (o^~  dk  a-ck,  bk) 
and  (a**i-r  6**,).  In  this  way  the 

wing  parts  arc  divided  into  wing  panels. 

It  is  convenient  to  have  the  panels  num¬ 
bered.  The  following  system  may  be  used. 
Panel  number  1  is  the  panel  in  the  first 
wing  part's  first  column  and  first  row. 
Number  2  stands  for  the  pane!  in  the  first 
column  and  second  row,  and  so  on.  r 
becomes  the  number  of  the  panel  in  the 
last  row.  The  first  panel  in  the  second 
column  has  number  r-fl.  The  last  panel 
of  the  first  wing  part  is  qr.  The  numbers  Tor 
the  other  wing  parts  follow  after  qr.  The 
number  i  of  a  panel  is  obtained  from  its 
location  in  the  wing  part  k,  panel  column  /. 
and  panel  row  m 

k-1 

i=mi(/-l)rtT  J  q,r • 

/-i 

where  q,  and  rt  arc  the  total  number  of 
columns  and  rows,  respectively,  for  wing 
part  j. 


Kg.  A 2.2.  Wing  pane)  with  horseshoe  vortex  and  con¬ 
trol  point. 

The  following  set  of  equations  can  be 
used  for  obtaining  the  quantities  that  specify 
the  vortices  xr,  yv.  It,  7.  and  the  control 
points  xD,  yD  from  the  quantities  that  deter¬ 
mine  the  panels  a,  b,  c,  d,  and  e.  Fig.  A  2.2 
illustrates  the  variables  for  vortices  and 
control  points. 

xv = Hr  ct  -  (3  dkm  B-kdK  u+1)/4  -f  (ijv~  bj)  -  tg  7. 
XD*=Qk-irC **04.rs"i‘3t4.  0*i)/4 

■V  [c*+l  ■  04*1.  a  3(4-,  1.  b-i)M 
—  Ck  ’  04.  n  4-  3t/i._  n- 1)/4  -f  ai+1  —  ak] 

*  (yv  bk)l(bk, ,  —  bj) 

b  —  O^.i-i  ~  c*.i)/2 

*S  ^  =  [cJfc*I  ‘  (3t4»i.  a  v  *4*1.  is-t-l)/d 

~  Ck(3dk.  a  T  4.  n*l)/l 

■V  °i+ 1  ~  J  ~  bk) 

APPENDIX  3 

THE  INFLUENCE  COEFFI¬ 
CIENTS  P,  Q,  R 

The  horseshoe  vortex  is  positioned  in  a 
plane  z  =  const,  with  the  bound  vortex 
forming  an  angle  to  the  y-axis  and  the 
trailing  vortices  parallel  to  the  .x-axis.  The 
geometry  is  obtained  from  the  real  wing  in 
a  Mach  number  M  air  stream  through  a 
Frandtl-Glauert  transformation,  that  is  the 
y  and  z  coordinates  arc  multiplied  by 

p,  p  =i/r:ip. 


Fig.  A2.1.  Wing  parts  and  panels. 
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The  three  disturbance  velocity  compo¬ 
nents  at  (lie  point  (xD,  pyD,  p:D)  due  to  a 
horseshoe  vortex  of  intensity  T  and  with 
the  geometry  defined  by  the  coordinates 
(xv,  ptjr,  ftzv)  and  the  dimensions  fib  and  ji 
will  be  calculated  by  Uiot-Savart's  law. 

The  bound  vorlex  forms  part  of  the  line 

iQpP'iy -Ur)  =  x-xv,  :  =  pzv. 

A  straight  line  vortex  induces  velocities  that 
lay  in  planes  normal  to  the  vorlex.  The  in¬ 
duced  velocity  in  a  point  P  is  also  normal 
to  the  straight  line,  which  connects  P  and 
the  vortex.  Fig.  A  3.2  shows  the  line  vortex 
between  Q  and  II  and  the  induced  velocity 
vector  in  P.  PS  is  the  line  through  P  normal 
to  the  vortex.  P'S  is  the  projection  of  that 
line  on  the  plane  z=pzv.  The  equation  for 
P'S  is 

~P  •  (y -yo)l  »g .«  «  x  -xD 


This  determines  the  point  of  intercept  5 

(>•  P*l>  P: r)- 

^  Xy-rXg  tg zp  +P-(yD- IJr)  tg/z 
’  tg^-r-i 

tsV^1 


Pn~- 


\Ya  will  denote  the  velocity  induced  in 


\V  _r(C0S/,~FC0S  y) 

4.i i  ¥+r- 


QS 

COS  y  =  =  = 

•  ah 


±  \  IP- C* 


QP  \A~Tlii^Cr~Zr 

cos  y  >  0  for  rj>  yr—Ii 
cos y  <  0  for  i)<yv—h 


ns 


COS  ff  =  ==  = 


± \' IP  -r  IP 


RP  I  Az^r  IP-IP-Z1 

cosy<0  for  ij>yv-rh 
cosip>0  for  j]<yv-rli 

a*- = (si>y  =  (xD-*y-  -  p'~(v  -yDy- 
B  =  p-(t]-yv  +  h) 

C  —  $  —xv  -i-  ph\gfl 

D  =  P  (tJi+lt-V) 

E  =  Xv-t  pil  tg  «  —  £ 

Z  =  p-(zy- :D) 

The  disturbance  velocity  1FB  is  resolved 
into  its  components  along  the  coordinate 
axis  uB,  vB,  wB 

uB  ——  IV^sinccos/t 
vB=  it  jsinesin/t 
iuB  =-\VBcose 
Z 


Slit  £  — 


cos  c  = 


+A 


VA--rZ- 


U-+z- 

cose>0  for  x0>£ 
cos  e  <  0  for  xD  <  £ 


The  port  free  vorlex  and  the  velocity  vec¬ 
tor  \VP  induced  it  P  arc  shown  in  Fig.  A 3.3. 


r  l+cos  - 


in  \l G-  +  Zz 


I;ig.  A  3.2.  The  bound  vortex. 


1G 
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The  conipoRcnts  or  Ws  along  the  axis  arc 
us.  Vs.  and  ws. 

us  =  0 

vs  =  Ujsint 
ws——  U’sCOST 


sinr= 


Z 


VP+Z* 


cosr= 


I 

\1 


COS  y.= 


H 


\,Gt-rH‘  +  Zs 


G=p-(yD-yr+K) 

i/=Xa-Xr-r/?htg/l 


The  components  of  velocities  along  the 
coordinate  axis  are  uF,  vr,  and  wF. 


tip  —  0 

vF  =—  Wrsinr 
Wp——\YpCosr 


sin  r* 


cos  »•= 


Z 


W±r- 

G 

l^Tz- 


The  velocities  Hrs  due  to  the  slarboanl 
free  vorlcx  are  obtained  similarly. 


rig.  A3.-I.  Starboard  free  vortex. 


T  cosa-f  1 

5st  fi*+z- 


coscr  = 


\i *Tj2+z2 


I  =  K'Jv  +  h~  Ud) 

J  =  XD  -  Xy  -  fill  tg« 


The  complete  horseshoe  vorlcx  causes  a 
disturbance  velocity  with  the  components 
u,  v,  w 

U  =  U  3  —  Up  -r  Us 
V  =  VB-T  Vp  —  vs 
W  —  tVg  T  Wp  T  IVg 


r  cosy -f  cosy  „ 

U— - /.  ;  ,■=;■  Z  COS  ll  =  P[' 

-lx  A--r/r 


y-rCOSns  .  1  -!-COS;< 

s'n«--^r— ^ 


r  fcosy-rco: 
cosa+n 

r-+r-  J 

r  r_cos-/-fcos<p  ,  1-fcosjj,.. 
W~M/~1?+ZS  Gr~r- 

_S£Z+>/].w 


(A3.1) 


The  rule  or  signs  for  the  first  term  in  in 
is:  use  minus  when  xD  >  f ,  use  plus  when 
~ld  >• 


APPENDIX  4 

EFFECT  OF  WING  ELASTICITY 

Eq.  (4)  yields  the  distribution  or  angle 
of  attack  that  corresponds  to  a  given  load¬ 
ing.  The  obtained  angles  for  the  loaded 
wing,  (oQ ,  arc  the  same,  whether  the  wing 
is  rigid  or  clastic.  However,  when  the  elastic 
wing  is  unloaded,  it  will  in  general  change 
its  shape.  The  corresponding  change  in 
slrcamwisc  angle  will  be  denoted  Aw.  It  is 
often  of  interest  to  know  the  distribution 
of  angle  of  attack  for  the  unloaded  elastic 
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wing,  «0.  Mostly  o>0  is  the  input  for  the 
calculation  of  the  loading  of  an  clastic  wing. 

{we}  =  -r  {A ca} 

The  angle  of  deformation  is  obtained 
from  a  deformation  matrix  [/)]  of  the  struc¬ 
ture  and  the  load  vector  {/}  for  the  de¬ 
formed  wing 

{Ao,}  =  [D]{1} 

The  load  vector  is  found  from  Eqs.  (5)  and 
(9)* 

=  [Sx]-'{o0} 

A  combination  of  these  equations  gives  the 


unloaded  wing  angle  of  attack  distribution 

w=fh  l&Hw 

LL  \J  **'  L  \J  J 

(A  4.1) 

The  loaded  wing  angle  of  attack  distribution 
is  obtained  fre’n  Eq.  (A  4.1) 

W 

(A  4.2) 

Inserting  Eq.  (A4.2)  into  Eq.  (5)  gives  an 
equation  for  calculation  of  loadings  for 
elastic  wings. 


